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The Fermi Motion and J/ip Production at the Hadron ColHders 

M.A. Gomshi Nobary, B. Nikoobakht 

Department of Physics, Faculty of Science, Razi University, Kermanshah, Iran. 

The Fermi motion effect in J /tj) production in various hadron colliders is studied. We deduce 
that in agreement with sum rules which the fragmentation functions should satisfy, while the effect 
has considerable influence on the fragmentation probabilities and the differential cross sections, the 
total cross sections are essentially left unchanged. 



I. INTRODUCTION 



II. FRAGMENTATION FUNCTIONS 
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The J /-ip state has been one of the interesting prob- 
lems of QCD both in theory [1,2,3] and in experiment 
[4]. QCD predictions for production cross section of 
J I'll} fails to agree with the experimental results. Color 
octet senario is introduced to bring about the agree- 
ment [5,6]. 

In this work, we introduce the Fermi motion into 
the J /ip production. We only include the charm frag- 
mentation contributions. The Tevatron Run I, the 
Tevatron Run II, the RHIC and the CERN LHC en- 
ergies are considered to demonstrate the effect. 



J 



With the choice of a light cone the wave function, 
we introduce the Fermi motion effect in fragmentation 
production of the J /ip bound state. In the leading or- 
der perturbative regime, the fragmentation functions 
for J /if production without and with the Fermi mo- 
tion are obtained [7]. 

In the absence of Fermi motion, where the confine- 
ment parameter is /3 = 0, the fragmentation function 
has the following form 



Dc^j/^{z,fio,(3 = 0) = 



z{l - zf Sf^z^ + 2^2^(4 - 4z + 5z2) 



16TOi^(z) 

-(16 - 32z + 24z2 - Sz^ 4- Qz^) 



(1) 



where as is the strong interaction coupling constant 
and Cf is the color factor. The quantity (fcr^) is the 
average transverse momentum squared of the initial 
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state heavy quark, the parameter ^ is defined as ^ 
(fcT^)/™^ and finally F{x) is given by 



F{z) = C'^^ -{z- 2)'(3z - 4) + Cz2(8 - 7z + z^) 



(2) 



r 



In the presence of Fermi motion the fragmentation function is obtained as 
I 



Dc^j/^{z,Ho,l3) 



TT^alClikT^y/^ f dqdx\i;M\''x^l - zfzq 



2m 



G{z) 



x^l-A{l-x)z + 2(4 - lOx + 7x^)z^ 

+A{-1 + - bx^ + 4x)z^ + (1 - 4.T + _ _^ ^i^^A 

+T]^z^ [l~2x + z^ + x^{2 - 2z + z^)] + ?7[2 + (-6 + 4:x)z 
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+(9 - 8a; + 2x'^)z'^ - 2(2 - a; + x^)z^ + (1 + x^)z^] 
+(,z^ [1 + 2a;^(2 - 2,z)z + z^ + Ix^z^ + 2a;(-l + - Iz^) 

W{2 - 80 + 9^^)] + rf'il - zf + f{l - xfx^z^y 

I 



(3) 



The function G{z) reads as 



{ [r?(l - zf + ix'z'' + (1 - (1 - x)zf\ + ^) + ^(-1 + x)xz'' - 1 + (1 - a; + x^)z\ }^ (4) 



D (z) xlO^ 



J/iA; /3=0 .2, 0.4,0.6 GeV 




Moreover here we have defined 77 = cf jm^ . 

The two cases of the Fermi motion off and on may 
be compared using the above fragmentation functions. 
Note that /3 = and 7^ correspond to these cases. 
For /3 value, we have chosen the range of /? =0.0 - 0.6 
GeV . The behavior of the fragmentation function (3) 
is shown in Fig. 1 for = 0, 0.2, 0.4 and 0.6 GeV. 
We take the maximum value of (3=0.6 GeV and use it 
for our further considerations. 



0.2 0.4 0.6 0. 



FIG. 1: Fragmentation function for J ftp production. 
While the dashed curve represents the function when the 
Fermi motion is off, the solid ones show the behavior of 
this function when the Fermi motion is on. The f3 values 

arc indicated. The function picks up as increases. This 
gives raise to increase in the fragmentation probability. 



III. INCLUSIVE PRODUCTION CROSS 
SECTION 



We have employed the idea of factorization to eval- 
uate the J/ip production cross section at hadron col- 
liders. For pp collisions we may write 



- — {pp ^ c^ J/iI;{pt)X) = ^ / dxidx2dzfi/p{xi,iJ,)fj/p{x2,lJ.) 

X [^{ij c{pt/z)X, ij,)D^^j/^{z,ij.,P) 

I 



(5) 



Where fi_j are parton distribution functions with 
momentum fractions of xi and X2, o is the charm 
quark hard production cross section and D{z,iJ.,(3) 
represents the fragmentation of the produced heavy 
quark into cc state with confinement parameter /3 at 
the scale fi. We have used the parameterization due 
to Martin-Roberts-Stiriling (MRS) [8] for parton dis- 
tribution functions and have included the heavy quark 
production cross section up to the order of [9] . The 



dependence on fi is estimated by choosing the trans- 
verse mass of the heavy quark as our central choice of 
scale defined by 



^iR = \/pT^ (parton) -|- 



(6) 



and vary it appropriate to the fragmentation scale of 
the particle state to be considered. This choice of 

scale, which is of the order of pr (parton), avoids the 
large logarithms in the process of the form ln(mQ / /u) 
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or ln(p'i7//). Howcvc!!', wc; have to sum up the loga- by the Altarelli-Parisi equation [10]. The following 
rithms of order of /Uij/wQ in the fragmentation func- form of this equation is used here 
tions. But this can be implemented by evolving them 

I 
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FIG. 2: The differential cross section for direct frag- 
mentation production of J/ip and its subsequent decay 
J/?/) M^M^ at the Tcvatron Run I energies. While 
the dashed curve is obtained using (1) or equally (3) with 
f3 = Q, the solid one is due to (3) with /3 — 0.6 GeV. 
The result is compared with the CDF Run I data. Other 
contributions are not included. The scale is chosen to be 



Here Pq^q{x = z/y,fi) is the Altarelli-Parisi split- 
ting function. The boundary condition on the evolu- 
tion equation (7) is the initial fragmentation function 
D^^jl^{z,ii,(i) at some scale = /Xo- In principle 
this function may be calculated perturbatively as a 
series in at this scale. 

Detection of final state requires kinematical cuts of 
the transverse momentum, pT, and the rapidity, y. 
We have imposed the required pSp'*" and y'^"'' in our 
simulations for different colliders as required and have 
used the following definition of rapidity y = \ \og[{E— 
Pl)/{E+pl)]. 



IV. RESULTS AND DISCUSSION 

In this work we have employed a light-cone wave 
function to introduce the Fermi motion in J /ip pro- 
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FIG. 3: The differential cross section for direct fragmen- 
tation production of J/ip at Tevatron Run II. The two 

curves arc obtained using (3) with /3 = and 0.6 GeV 
respectively. The scale has been set to 2/ii{. 



duction in direct fragmentation channel and obtained 
its fragmentation function in leading order perturba- 
tive regime. The behavior of J/f/' fragmentation func- 
tion with such a wave function is illustrated in Fig. 
1. 

We have studied the production rates at which the 
J /ijj state is produced with and without the effect 
of Fermi motion. First we present the pt distribu- 
tion of BR(J/'!/) — > fX^ ij,~)d(j / dpT for the cases of 
the Fermi motion off and on along with the CDF 
data at Run I in the Fig. 2. The branching ra- 
tio BR(J/'^i' /u+/i-)=0.0597 is taken from [11]. 
The poor agreement with data is due to the fact 
that here we have only considered the contribution 
oipp ^ J/tp ^ fJ,'^H~. Similar behavior at Teva- 
tron Run II energies is shown in Fig. 3. We have 
also extended our study to the cases of the RHIC and 
the CERN LHC pp colliders. Here we provide the 
Pt distributions of the differential cross sections for 
J/tp and compare the two cases of the Fermi mo- 
tion off and on in the Fig. 4. In all cases we have 
used /3 = 0.6 GeV for the confinement parameter in 
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FIG. 4: The differential cross section for direct fragmenta- 
tion production of J/ip at the RHIC and the CERN LHC. 
The shift due to the Fermi motion is significantly increased 
in the case of LHC but shows to be less important at the 
RHIC. In both cases the fragmentation function (3) is em- 
ployed with /3 = and 0.6 GeV respectively. 



the fragmentation functions. Naturally, the results for 
(3 in the range of - 0.6 GeV fall between the above 
results. 

Calculation of the total integrated cross sections 
for each case shovirs that the total cross sections for 
with and without Fermi motion essentially remain un- 
changed within the uncertainties of Mote Carlo sim- 



ulations. The reason is first due to the momentum 
sum rule which the fragmentation functions should 
satisfy. In other words although the modification of 
fragmentation functions by the Fermi motion redis- 
tributes the final states, the integrated cross sections 
are left unchanged. Alternatively although the Fermi 
motion increases the fragmentation probability for the 
state, i.e., introduces a state with overall higher mass, 
the cross section is lowered by just the same amount 
when we introduce the effect in calculation of the total 
integrated cross section. It is evident from the Figures 
2,3 and 4 that the effect increases with increasing ^/s. 
The kinematical cuts play important role apart from 
V^. The large cross section at the RHIC compared 
with the LHC in the Fig. 4 is an example. 

There are two main sources of uncertainities.The 
first is about the simulation of J /if) production at 
hadron colliders such as the uncertainties along with 
the fragmentation functions and parton distribution 
functions. These kind of uncertainties are well dis- 
cussed in the literature. The second source of uncer- 
tainty is due to the choice of the confinement param- 
eter in the fragmentation function. Relying on our 
discussion in section 2, our choice of /3 = 0.6 GeV 
seems to be justified. Future determination of this 
parameter will shed more light on this situation. It 
is worth mentioning that our choice of charm quark 
mass, i.e., 1.25 GeV, have put our results in their up- 
per side and that the change of the charm quark mass 
in its acceptable range does not have significant im- 
pact on the Fermi motion effect in J /tp production. 

More detailed information concerning this work ap- 
pears in [12]. 
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